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ABSTRACT 

We present a quantitative model of the magnetic energy stored and then re- 
leased through magnetic reconnection for a flare on 26 Feb 2004. This flare, well 
observed by RHESSI and TRACE, shows evidence of non-thermal electrons only 
for a brief, early phase. Throughout the main period of energy release there is a 
super-hot (T ^ 30 MK) plasma emitting thermal bremsstrahlung atop the flare 
loops. Our model describes the heating and compression of such a source by lo- 
calized, transient magnetic reconnection. It is a three-dimensional generalization 
of the Petschek model whereby Alfven-speed retraction following reconnection 
drives supersonic inflows parallel to the field lines, which form shocks heating, 
compressing, and confining a loop-top plasma plug. The confining inflows pro- 
vide longer life than a freely-expanding or conduct ively- cooling plasma of similar 
size and temperature. Superposition of successive transient episodes of localized 
reconnection across a current sheet produces an apparently persistent, localized 
source of high-temperature emission. The temperature of the source decreases 
smoothly on a time scale consistent with observations, far longer than the cool- 
ing time of a single plug. Built from a disordered collection of small plugs, the 
source need not have the coherent jet-like structure predicted by steady-state re- 
connection models. This new model predicts temperatures and emission measure 
consistent with the observations of 26 Feb 2004. Furthermore, the total energy 
released by the flare is found to be roughly consistent with that predicted by 
the model. Only a small fraction of the energy released appears in the super-hot 
source at any one time, but roughly a quarter of the flare energy is thermalized 
by the reconnection shocks over the course of the flare. All energy is presumed 
to ultimately appear in the lower-temperature T ^ 20 MK, post-flare loops. The 
number, size and early appearance of these loops in TRACE'S 171 A band is 
consistent with the type of transient reconnection assumed in the model. 
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Introduction 



Si nce the introduc tion of magnetic reconnection iGiovanellil ( 119471 ); IS wee tl (119581 ); iParker 



( 119571 ) ; |Petschekl (ll964[ ). many model s have propo s ed it as an eleme n t of fl ar e- related phenom - 



ena, such as chromo s pheric ribbons ICarmichaell ( 119641 ); ISturrock (196 81); iHiravam, 



Kopp and Pneumanl (119761) . associated coronal mass ejections 



Moore and Labont 



Forbes and Priestl ( 19951 ). above-th e-loop X-ray sources Masuda et al. ( 1994 ); 




(1997) and supra-arcade downflows lMcKenzie and Hudson! (119991 ) ; iMcKenzid (120001 ) ; ISheeley. Warren, and 
(120041 ). 



The term "magnetic reconnection" is commonly used to refer to either or both of two 
different effects: topological changes to magnetic field lines and conversion of magnetic energy 
into other forms. Only the first effect, topological change, requires a parallel electric field and 
thereby qualifies as "magnetic reconnection" in the strictest sense. It has been successfully 
quantified by several kinds of observation. The amount of flux being top ologically changed 



can be computed by tracking the ri bbons acro s s phot ospheric field in a flare Forbes and Priest 



(I1984T): IPoletto and KoppI (119861): biu et all (120021 ). computing the flux ejected by CMEs 
Qiu et al\ (120071): lLongcope et all (120071 ). or measuring the fluxes of retracting downflows 
McKenzie and Savage! (120091 ) . Energy conversion, on the other hand, can occ ur together wit h 
the topologi cal change, as in the diffusion-dominated Sweet-Parker models IParker! (119571 ): 
Sweetl (119581 ) or as a subsequent response to it, as in the Petschek model. The latter case 
would be more properly termed "post-reconnection energy conversion" , but Petschek's model 
is customarily referred to as one of reconnection. 

Recent theoretical progress has revealed that in order to be genuinely fast, the re c onnec - 
tion electric field must be lo c alized to a smal l portion of a current sheet iBirn et all (120011 ): 
Biskamp and Schwara (120011 ): iKulsrudl (120011 ). Since such a small region will contain negli- 
gible magnetic energy, energy conversion must occur away from the region of flux transfer, 
as a post-reconnection response to small-scale topological change. This requirement of phys- 
ical separation between the non-ideal electric field and the energy conversion offers the hope 
that energy conversion can be modeled without direct knowledge of, or appeal to, the specific 
physics dictating flux transfer. 

In Petschek's prototypical model of fast magnetic reconnection, energy conversi on occurs 



(1964 


); 


Sonnerup 


( 


1970 


); 


Vasvliunas 


( 


1975): 



els invoked a resistive electric field to transfer the flux, it has since been established that 
any localized process will produce the same external response provided the plasma may 
be treated by fluid equations (and the reby support collisional shocks) on its largest scales 



Erkaev. Semenov. and Jamitskyi (120001 ). The two-dimensional, steady-state model has also 
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been 



(2001 



1982 



generalized to transient reconnection iBiernat. Hevn. and Semenovl (I1987D; iNitta et al. 



, reconnection between skewed mag netic fields (jPetschek and Thornei 11967 



Soward 



Skender. Vrsnak. and Martinis! 12003k i.e. with a magnetic field component, or "guide 



field" , in the ignorable direction, sometimes c alled "two-and-a-half" dimension s ) and f inally 
to transient reconnection in three dimensions lLongcope. Guidoni. and Linton! (j2009bl ). We 
find that all such models predict very similar energetics since they all model external fluid 
responses to an assumed flux transfer of negligible dimension. 



Tsuneta (11996! ) presented a thorough application of Petschek's reconnection model to 
observations of a solar flare. He demonstrated that the temperature structure observed above 
the long-duration event on 21 Feb 1992 was consistent with the SMSs of two-dimensional 
Petschek reconnection. He showed further that local energy fluxes were consistent with 
the model, but could not, for a number of reasons, turn this into a global accounting of 
net energy released. First of all, the flare took place on the eastern solar limb making a 
reliable model of the coronal magnetic field difficult. Indeed, invoking Petschek's original 
two-dimensional model tacitly neglects any possible magnetic-field component along the line 
of sight. Secondly, the reconnecting current sheet was probably a dynamical after-effect of 
a CME whose energetic contribution could overwhelm that of the flare itself. Finally, the 
impulsive phases of most large flares appear to include non-thermal particles, which fluid- 
based reconnection models ca n not easily incorpor a te, but which accou nt for a significant 
energy iLin and Hudson! (1197ll ): Istrong et all (119841 ); Emslie et all ( 120041 ) . 

Except for the adverse geometry, the difficulties faced by Tsuneta are likely to 

be common to most large flares. Non-thermal particles pose an especially difficult challenge 
for global energetics. Wh ile much progress has been made und erstanding the mechanism 
for the their acceleration ( jMiller et all 119971 ; lAschwandenl |2002| . for reviews), most models 
pre-suppose a non-ideal electric field, fluid turbulence, or plasma waves. Large-scale fluid 
models of reconnection are consistent, to some extent, with each of these, but do not include 
coupling to non-thermal particles. Since these particles ultimately account for a significant 
part of the released energy their effect is almost certain to be non-negligible. 

A small subset of large flares show little evidence for non-thermal particles, but instead 
exhibit hard X-ray (HXR) thermal bremsstrahlung from a super-hot plasma (T ^ 30 MK). 
Rare as they are, these flares offer a unique opportunity to comprehensively quantify the 
energy release in a flare using a fluid reconnection model. Furthermore, these flares are 
often of the so-called compact variety and thus lack CMEs which could otherwise muddy 
the energy estimate. Lacking an eruption, the pre-flare magnetic field can be assumed to be 
near equilibrium and its energy more easily estimated. Such an ideal flare occurred on 26 
Feb 2004 and was well observed by both RHESSI and TRACE at high cadence with little 



-4- 



interruption. This flare has been found to be one of the rare (one in ten) X-class flares 
unassociated with any CME IWang and Zhang (120071 ). We present here a model of how 
magnetic reconnection produced the phenomena observed in this flare, including estimates 
of the entire energy release process. 

Spectrally resolved HXR images show a compact source probably located atop the flar- 
ing loops. This source is assigned an electron density which is high but typical of flaring 
plasma, (n e ^ 10 11 cm -3 ). The high density is required by both the large observed X-ray 
flux and in order that electrons of such high energies to be collisional y thermalized . High - 
temperature loop-top sources such a s this h ave been ob s erved before lActon et all (119921 ); 
Petrosian. Donaghy. and McTiernanl (120021 ); iJiang et al\ (120061 ). and have proven particu- 
larly challenging to explain. We believe that shocks from fast reconnection offer a promising 
avenue for modeling the source, since they would both compress and heat the plasma. 

The source's persistence time vastly exceeds the time that it would take for the localized 
structure to conductively cool or to freely expand. Some models have addressed the former 
difficulty by proposing that a turbulent reduction of the electron thermal c onductivity per- 
mits the high-temperature plasma to cool more slowly IJiang et al\ (120061 ). This does not 
explain, however, how its high pressures are confined to the top of the loop. A second ap- 
proach is to assume the observed, persistent structure is actually a super- position of s equen - 
tial transient events ; this a pproach is common in models of flar i ng loops iHori et al. ( 19971 ) ; 
Reeves and Warrenl ( 20021 ); IWarren. Winebarger. and Hamilton! ( 20021 ) ; IWarrenl ( 20061 ). The 
temperature of these features does, however, change smoothly and gradually, in a manner 
resembling a slow cooling. Proposing that it consists of multiple isolated elements then re- 
quires an explanation for their smooth collective behavior. This is presumably a signature 
of the evolving energy conversion process — possibly magnetic reconnection. 

In this work we show that magnetic reconnection occurring sporadically in small patches 
across a pre-flare current sheet will produce a super-position of high-density, high-temperature 
sources consistent with the observations on 26 Feb 2004. Moreover, the values of the tem- 
perature and emission measure, as well as their gradual temporal evolution, is consistent 
with the model. The specific physics responsible for transferring flux within the reconnec- 
tion patches does not enter into these predictions and cannot therefore be constrained by 
the observations. The mean rate of transfer, which is to say the number of new patches per 



unit t i me, does enter, but this is a global quantity that can be measured iForbes and Priest 



(j 19841 ); IPoletto and Koppl (119861 ); iQiu et all (120021 ). While this kind of disordered process is 
less obviously related to Petschek's model than a long, steady outflow jet, we show that it 
is consistent with the observed size and temperature of the loop-top source. 



This combination of observations and model are presented as follows. The next section 
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describes the observational data of the flare. Section 3 then presents a model of the pre- 
flare magnetic field, including the geometry of the three-dimensional current sheet on which 
magnetic reconnection subsequently occurs. We then present, in section 4, a general model 
of transient, localized magnetic reconnection within a current sheet. The Alfven speed 
and angular discontinuity characterizing the pre-reconnection current sheet directly predict 
the temperature of the reconnection outflow. We show that this prediction is common 
across all models of fast reconnection: steady state or transient, two-and-a-half or three- 
dimensional. In section 5 we apply this model to the observations of 26 Feb 2004. We 
maintain that evidence in this particular flare favors a transient reconnection model. We 
predict an emission measure that can be directly compared to observation. We use different 
means to measure the mean flux-transfer rate and find that all agree and predict flaring 
emission measure comparable to the observation. Section 6 uses the reconnection to model 
the full flare, including an account of its energetics. Finally, section 7 discusses the possible 
applicability of the new model to a broader class of flares. 



2. The Flare: 26 Feb 2004 

2.1. Active region 10564 — emergence 

The 26 Feb 2004 flare (SOL2004-02-26T02:03:00L161C076) occurred within the com- 
plicated, but basically bipolar, active region 10564. The simple bipolar active region had 
one umbra of each polarity until, on 24 Feb, new flux began to emerge into it. This new 
flux eventually developed its own umbrae and the new positive and old negative penumbrae 
merged (see Figured]). After roughly 50 hours of steady emergence an X-class flare occurred 
within the region. 

To distinguish between new and old flux we partition the line of sight MDI mag- 
netogram as shown in Figure [2J The radial component, B r , is derived from the line- 
of-sight measurement assuming the field to be perfectly radial^] Those pixels exceeding 
\B r \ > B t hr = 125 G are grouped into uni p olar regions according to the horizontal gr adient 
in B r iBarnes. Longcope. and Lekal (120051 ); lLongcope. Barnes, and Beveridgd (j2009al ). The 



regions labeled N01 and P03 constitute the old bipole and all other partitions are from 
emergence since the beginning of 24 Feb. 



1 To assess the accuracy of the radial-field assumption Figure[T]has a longitude axis and marks showing the 
time the old polarities cross central meridian. Any systematic error would lead to an artificial center-to-limb 
variation centered at the x . 
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Fig. 1. — MDI data of new flux emergence in AR 10564. Rows are data from four successive 
times, progressing upward, spanning 54 hours. Printed on each panel is the time since 
00:00UT on 26 Feb 2004. The left column is the line-of-sight magnetogram scaled between 
±1000 G; the right is the continuum image. 
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Fig. 2. — Partition of the line-of-sight magnetogram from just before the flare. The lower 
panel shows the magnetogram, scaled to ±1000 G. Lines denote the boundaries of a partition, 
and + and x are the centroids of each region. The upper panel shows the same centroids 
along with null points (triangles) and the principal spines (solid) connecting them. The 
dashed curves are separators connecting null point Al to B3 and B4. 
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The net flux and centroid of each region ($j and Xj) is found by integrating its surface 
field. The time evolutions of old flux, $po3 (positive) and &noi (negative), are shown as 
dashed curves in Figure At the beginning of the data each has approximately 8 x 10 21 Mx 
of flux. Over the course of our data set the fluxes in each polarity gradually decreases until 
they are $ P0 3 = 5.6 x 10 21 Mx and § N01 = 7.4 x 10 21 Mx at the time of the flare. The new 
flux, comprising the remainder of the regions in the partition, increases much more rapidly 
over this period reaching twice that of the old flux (|$| ~ 13 x 10 21 Mx) by the time of the 
flare. The emergence is steady and there is no clear feature in the evolution at the time of 
the flare (vertical dashed line). There is a small excess of negative flux (broken line) due to 
the omission of positive flux either outside the field of view or below the threshold strength, 
.Bthr = 125 G. 
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Fig. 3. — The fluxes in the new (solid) and old (dashed) polarities of AR 10564; the broken 
line is the total signed flux of all regions. The value at the time of the flare, 26 Feb 2:03UT 
(vertical dashed line), is given and indicated with a symbol. The signed sum of all fluxes, 
shown as a broken line, is close to zero. The times at which each of the old polarities crosses 
central meridian is shown with an x, and the longitude of their mid-point is shown by the 
bottom axis. The top panel shows the GOES lA - 8A flux. The top axis gives hours past 
24 Feb 00:00. 

We interpret flare observations in terms of a magnetic model built using these photo- 
spheric observations combined with assumptions about flux emergence. The partitioning of 
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the photospheric field described above makes no assumption about the state of the coronal 
f ield, but does provide a framework for quantifying its connectivity lLongcope and Klapper 
( 120021 ). Regardless of the state of the coronal field, every coronal field line can be assigned to 
a domain based on the photospheric regions of its footpoints: its connectivity. Different coro- 
nal fields have different connectivity and the lev el of difference between two possi b le field s 
is quantified by the differences in domain fluxes lLongcope. Barnes, and Beveridgei (j2009al ) . 
We use the potential field extrapolated from the pre-flare magnetogram of 1:39 (Figure [2]) as 
a reference for constructing such differences. The potential field has the minimum possible 
magnetic energy so the d egree by which a given field differs from i t relat es to the amount 
for free energy in the field lLongcope] (120011 ) ; lLongcope and Magaral (120041 ) . 



Provided the post-flare coronal field has lower magnetic energy than the pre-flare field, 
we expect its connectivity to be closer to that of the potential field. The potential field 
conne ct ivity is characterized in term s of the magnetic skeleton iPriest. Bungey. and Titov 
( 1997 ); Longcope and Klapper ( 2002 ); Beveridge and Longcope (j2005l ). illustrated by the 
top panel of Figure [2j This stylized figure shows sources (+s and xs) and several of the 
magnetic null points (triangles) of the potential field. Solid curves show the spine lines from 
each null point, and dashed lines are two separators which extend into the corona from null 
points Al, B3 and BA. The spine lines trace out the edges of separatrix surfaces which divide 
the coronal magnetic field into different domains. Most significant for the present case, they 
divide the pre-existing coronal field, anchored to P03 and N01, from newly emerging field. 

The pre-flare field will differ significantly from the potential field since it must contain 
significant currents. Rather than try to measure these currents we infer them from the con- 
nectivity after assuming that no significant reconnection occurred prior to the flare. This 
means that in the pre-flare field the old photospheric regions, P03 and N01, remain uncon- 
nected to the others, such as P09 or iV07, which recently emerged. The potential field, and 
therefore the post-flare field, contains significant connections between new and old regions 
regions P09 and iVOl, due to their proximity. By our flux-emergence assumption, however, 
this connection contains zero flux in the pre-flare field. It is differences such as this that lead 
to large free energies in the pre-flare field. The forging of new connections through recon- 
nection will reduce this difference thereby releasing stored magnetic energy and presumably 
powering the flare. 



2.2. The flare loops 



EUV images provide evidence supporting our assumption that flare-related reconnec- 
tion creates field lines with new connectivity. TRACE obtained 171 A images at a 30-second 
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cadence throughout the flare. FigureHJshows an image from before the flare (1:31) and seven 
subsequent images during the flare. Post-flare loops occupy the region between by the spine 
curves, corresponding to a coronal volume beneath the associated separatrices. In images 
from early in the flare the 17lA emission is primarily chromospheri c, including moss, foot- 
points and flare ribbons (e.g. 1:55:14). Following a common pattern iGorbachev and Somov 
( 119881 ) : lLongcope et al\ (120071 ). the flare ribbons follow the topological spines. Coronal loops 
are visible as early as 1:57 and become the primary features by 2:04:14. 




Fig. 4. — TRACE 171A images from times before (1:31) and during the flare, shown in 
inverse greyscale. The grey scale is enhanced by a factor of eight in the first image and 
by 4/3 in the final three (2:14 - 3:01). Dark lines in each frame show the spines from the 
magnetic skeleton from Figure |2j 

The visible coronal loops are summarized in Figure [5j 143 loops were visually identified 
and manually traced in the collected 30-second cadence 171 A TRACE images. Those loops 
persisting longer than the median lifetime, 6.0 minutes, are plotted over the partitioned 
magnetogram with color indicating the time of first appearance. Early, short loops connect 
new positive regions P09 and P04 to the old nearby region A r 01. Later in the flare, longer 
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connections are made between iVOl and the more distant region P02. This is still newer 
than region P03 which was part of the original bipole. The apparent lengths of these loops 
increase with time approximately linearly according to the empirical relation 

LeuvW = L(t-t ) + L , (1) 

where L = 4.7 Mm hour -1 and L$ = 6.3 Mm using to =2:00:00. At the onset of the flare 
new loops were appearing at a rate of four per minute, as indicated by the dotted line in the 
bottom panel. 

The overall time history of the flare is summarized in the composite light curve of Figure 
|6j The first evidence of the flare is an impulsive brightening, at 1:53:50, in high energy X- 
rays (RHESSI 25 - 50 keV channel at top). The lower energy X-ray channels begin rising 
at about this same time, as do both of the GOES X-ray channels. The first 171A image 
after a five-minute gap (1:55:14) has total emission clearly elevated from the pre-flare level. 
Three minutes after the initial burst, at 1:56:40, all X-ray channels rise more rapidly to a 
higher level. It is during the times following this second rise that TRACE 17lA images are 
dominated by loops. 



2.3. Super-hot thermal loop-top source. 

Hard X-ray emission from the flare occurs in two phases with distinct spectral forms 
typified by the examples shown in Figure [71 Each spectrum is a 20-second integration of 
RHESSI detector 3. Plasma properties such as tempera ture and emissio n measure are found 



by fitting the spectra with the OSPEX software package ISchwartzl ( 119961 ) . In determining the 
best-fit model for each spectrum examined here, the following considerations were taken into 
account: spectral analysis from a single detector at a time, individual background subtraction 
for each energy range (6 - 12, 12 - 25, 25 - 50, 50 - 100 keV), albedo correction, the presence 
of instrumental lines at 8.5 and 10.5 keV (not shown in the figure), modifications necessary 
due to individual detector response differences and pulse pile-up and the attenuator change 
from Al (thin only) to A3 (thick and thin) at ^1:57 and back to Al at ^2:07. 

For every 20 sec. time interval from 01:54 - 02:10, we modeled spectra in two distinct 
ways. First we fit a thermal bremsstrahlung plus a non-thermal thick target (power law). 
Next we did a separate fit to two thermal bremsstrahlung components (a two-temperature 
fit). We found that the first phase of the flare (approximately 01:53:50 - 01:56:40) is best 
characterized by the first fit, a thermal plus non-thermal spectral model; we refer to this 
phase of the flare as the non-thermal phase. The second phase (approximately 01:56:40 - 
02:10) is best fit by the two temperature model with no non-thermal contribution; we refer 



2:00 2:15 2:30 2:45 3:00 3:15 

Fig. 5. — Coronal loops visible in TRACE 171A images are shown as color curves atop 
the partitioned MDI magnetogram. Colors from red to green represent times form 2:00 to 
3:00UT, as indicated in the top color bar. The upper panel shows the loops which appear 
before 2:18, and the middle panel shows later loops. For clarity, only those loops visible for 
longer than the median lifetime are shown. The bottom panel show the accumulated number 
of loops visible before any time. Horizontal bars show the lifetime of each loop, colored with 
the same code as the loops themselves. 
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Fig. 6. — A composite light curve of the flare on 26 Feb 2004. All curves show intensity 
on a linear scale. The bottom curve is the integrated intensity of TRACE 171A from the 
region shown in Figure HI The zero level is suppressed, but the level prior to 1:50 shows 
the pre-flare background. The dashed section is a 5-minute data gap and diamonds show 
the times at which the central six panels of Figure HI Above this are the curves from the 
low energy (Lo, 1 - 8A) and high energy (Hi, 0.5 - 4A) channels of GOES; the latter has 
been multiplied by two to facilitate comparison. The dashed c urve above th em is the plasma 
temperature derived from the ratio of GOES energy channels iGarcial (11994 ). Top curves are 
integrated intensity from RHESSI channels, from bottom to top, 6 - 12, 12 - 25 and 12 - 25 
keV. Vertical dotted lines show the beginning and end of the non-thermal phase of the flare. 
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to this phase as the thermal phase. During the thermal phase, the higher energy end of 
the spectrum (T ^ 12 keV) is consistently better fit by a super-hot thermal bremsstrahlung 
component than by a thick target power law component. In the 01:58:10 spectrum in Figure 
[TJ the super hot component (magneta) has a temperature of 46 MK 




Fig. 7. — RHESSI spectra from detector 3 integrated over 20 seconds plotted in black and 
multi-component fits. Separate components of each fit are plotted in different colors, and 
their sum is the red curve. The fit was performed over the region delineated by the vertical 
dashed lines, (left) An early phase in the flare, centered at 1:54:50, fit by one thermal 
component (T = 22 MK, blue) and thick-target emission from a non-thermal component 
with power-law spectral index 5 = 6.4 (green), (right) A later phase in the flare centered at 
1:58:10 fit by two thermal components:T = 46 MK (magneta) and T = 23 MK (blue). 

The spatial morphology of X-ray emission also differs between the two phases, as illus- 
trated by the RHESSI CLEAN images shown in Figure [8J During the non-thermal phase 
(Figure Ek) the highest energy emission (25 - 50 keV, green) originates in two distinct sources 
located on opposite sides of the PIL. These are presumably from non-thermal particles strik- 
ing the chromospheric footpoints of the reconnected field lines, perhaps linking P09 - iV06. 
While the lower energy emission (6 - 25 keV, red and blue) extends over the footpoints, it 
is more concentrated between the two sources. This point is presumably closer to the apex 
of the loop, or the site of the reconnection. It is notable that the peak of the 6-12 keV 
emission (red) occurs near the separator field line (magenta). 

During the thermal phase (Figures [Sb - d) the high and low energy emission originate 
from a single source near the separator field line. It is elongated approximately parallel to 
the separator, with length « 16 Mm and width « 9 Mm. Based on the magnetic geometry 
we believe this source is at the top of the post-flare loop system; it is a loop-top source 
viewed from above. By 2:05 (Figure [HH) the count rate in the highest energy channel (25 - 
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Fig. 8. — RHESSI images formed using the CLEAN algorithm on 20-sec integrations of 
multiple detectors. The grey scale in all images is from low-energy channel (6 - 12 keV) 
integrated beginning at 1:54:00. Contours are at levels of 50%, 70% and 90% of maximum 
for energy levels 6-12 keV (red), 12 - 25 keV (blue) and 25 - 50 keV (green). The panels 
show integrations beginning at A: 1:54:40, B: 1:58:00, C: 2:02:00 and D: 2:05:00. Also shown 
are the spines (black), null points (triangles), one separator (magenta) and polarity inversion 
line (violet dashed). 
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50 keV) had fallen too low to permit accurate imaging. 

Spectra were formed from both detectors 3 and 4 over each 20-second interval after 
1:54:40. All spectra resemble one of the typical cases shown in Figure [TJ but with different 
fitting parameters. Time histories of the best-fit temperatures and emission measures of 
each component are shown in Figure [9j Spectra at times before 1:56:40 have non-thermal 
(NT) components and are best fit by a thick target power law and one thermal component 
whose temperature and emission measure are plotted alone. All later spectra are fit by two 
temperatures. By 2:07:00 the two temperatures are sufficiently close that it becomes difficult 
to discriminate between the components. As a result, the fits from the different detectors 
give notably different values as indicated by larger error bars. As early as 2:05:00, the time 
of image EH, the temperature of super-hot component had fallen sufficiently low (26 MK) 
that it does not contribute enough signal to the 25 - 50 keV channel to permit accurate 
imaging; for this reason there is no green contour in that image. 

The super-hot component reaches as high as T = 48 MK at the beginning. It decreases 
on an approximately exponential decay (dashed curve) with a time constant of 14.2 minutes, 
while the cool component decays even more slowly. It is notable that the cool component 
matches the properties derived from the two GOES channels (after subtracting a background 
level determined from an interval prior to 1:00); the GQES-de rived emission measu re appears 



to be larger than the RHESSI value by roughly a factor of 2 lHannah et all ( 120081 ). 



3. Model of the pre-flare magnetic field 

The picture which emerges from the foregoing data is of a sudden burst of reconnection 
between the old bipole and the newly emerg ed flux. In the seminal model of such recon- 



nection, proposed by Heyvaerts et al. (119771 ). a current sheet forms between the two flux 



systems. Reconnection across this sheet releases the energy which powers the flare. A po- 
tential field above that configuration would contain a single coronal null point from which 
four separatrices divide the coronal flux into four domains: new flux, old flux and two mixed 
old/new domains. Remaining potential during emergence requires reconnection across the 
null point in order to convert new and old flux into the mixed type field lines. Without 
reconnection a current sheet forms at the erstwhile null point. 



The Heyvaerts et al. (119771 ) model is two-dimensional and the emerging field is purely 
anti-parallel to the overlying field. In the present case, however, the emerging flux is largely 
parallel to the existing bipole P03 - iVOl. This leads to a modified three-dimensional model 
in which the potential field contains a separator rather than a coronal null point. A potential 
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Fig. 9. — The temperatures and emission measures of the two thermal components from 
detectors 3 and 4. Curves show the simple average of temperatures from the two detectors 
(bottom) and geometric means of emission measures (top). Error bars show the values 
from each detector. The smoother solid curves are the values derived from the two GOES 
channels. The dashe d cur ve is the GOES-derived emission measure divided by a factor or 2, 
after Hannah et al. (120081 ). Lighter points show times after 2:07, where the two components 
were more difficult to discriminate. The dashed temperature curve is an exponential fit to 
the super-hot component; its decay time is 14.2 minutes. The imaging intervals, A - D, from 
Figure |8] are indicated along the bottom. 
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field extrapolated from the magnetogram of Figure [2] contains no coronal nulls but many 
separators; two separators are shown in the upper panel. Each separator is a magnetic 
field line (B 7^ 0) lying at the intersection of two separatrices. Together these separatrices 
separate four different field line connectivities and the separator lies between all four. 



In analogy to the two-dimensional model of Heyvaerts et al. (119771 ). emergence without 
reconnection leads to a current ribbon approximately following the path of the potential-field 
separator. Figure [TU] shows the configuration of such a current ribbon along the separator 
connecting photospheric null points Al and B4. There is a component of magnetic field 
parallel to the current, often called a guide field. In spite of its three-dimensional structure the 
equilibrium stru cture of this current rib bon ca n be approximated using the two-dimensional 
model of Green (jl965[ ) and Syrovatskii (Il97l[ ). 



Fig. 10. — A perspective view of a current ribbon along the separator field line connecting 
null points B4 and Al (magenta). The current sheet is indicated by green ribs and the 
current flows roughly parallel to the magenta separator. Examples of post reconnection field 
lines are indicated by red (P04 - N01) and blue (P09 - iV05). A blue asterisk shows a point 
along the separator located z = 7.9 Mm above the photosphere. 

To analyze the current ribbon, consider a plane perpendicular to the separator field line 
at one single point, such as the blue asterisk in Figure [101 Take this as the origin of the plane 
with local coordinates y and z. In the neighborhood of the origin the potential magnetic 
field within this plane is 

B^{y,z) = B g x - B' ± (yz + zy) , (2) 



where B' ± is the local field gradient at the origin lLongcope and Silval (119981 ) ; lLongcope and Magara 
J2004T ). At the coronal point indicated in Figure H0| B g = 331 G and B' ± = 34GMm _1 . 



- 19 - 



The plane's coordinates have been oriented so the field is radial along diagonals 45° 
from the coordinate axes. Unlike the two-dimensional case, these diagonal lines cannot be 
associated with the field's separatrices, since those are global features not derivable from 
local properties of the field. 

Th e perp endicular extent of th e current sheet follows the two-dimensional model of 
Greene ( 119651 ) and Syrovatskii (jl97ll ). The net current in the sheet, /, determines its width 



A 




(3) 



where B' ± > but / may be of either sign. The sheet extends along the ±y direction when 
I > and along ±z when / < 0, as shown in Figure HU Since B' ± varies along the separator 
so does the width of the current ribbon, as evident in Figure [TDJ It does, however, carry the 
same current / over its entire length, L. 



end-on view 



side view 




Fig. 11. — The geometry of a current sheet with a guide field component. The left panel 
shows the conventional "end-on view" of the sheet, along its current — x. Two Y-type neutral 
points, separated by A define the ends of the sheet. Viewing this from the left provides the 
less- conventional "side view", shown in the right panel. The Y-type neutral lines extend 
horizontally along the top and bottom of the sheet (dotted lines). The guide field is directed 
to the right (+x) on both sides of the sheet. The field on the near side (solid) is directed 
upward while that on the far side (dashed) is directed downward. The angle between these 
field lines is A6. 

The current sheet arises from a hypothesized absence, during pre-flare build-up, of 
reconnection required to achieve a potential field state. A potential field extrapolated from 
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the magnetogram of Figure [21 contains field lines linking, for example, the new pole P04 to 
the old pole N01. In fact roughly 60% of the flux from P04 connects to iVOl for a total flux 
^4-! = 2.8 x 10 21 Mx, where the superscript is a reminder that this value pertains to the 
potential (vacuum) field. 

All new polarities presumably emerged connected only among one another. In this case 
a connection between new and old regions would contain zero flux, \&4-i = 0, in contrast to 
the non-zero flux in the potential field. It is this discrepancy 

A* 4 -i = *4-i - , (4) 

between the p otential fiel d conn ectivity and the actual connectivity which leads to the sep- 
arator current iLongcopd (J200l|). For small values, the net current follows a simple self- 
inductance relation 

J A* 

c ~ SI ' < 5 » 

where we have omitted, for clarity, a factor depending on ln(|/|). 

The sense of current depends on the sense of change, A\I/, relative to the sense of the 
potential field. Potential field links the separator of Figure [10] in a positive sense: right- 
handed following the separator in the direction of the guide field. Before any reconnection, 
^4-i = 0, so the self flux from the current ribbon must be directed through the separator in 
the opposite sense. Thus the current must flow anti-parallel to the guide field: / < 0. This 
produces the roughly vertical sheet shown in Figure [101 

The separator shown in Figure UM has length L = 54 Mm and will therefore, according 
to eq. fl5J), carry I/c ~ 4 x 10 10 Gcm (i.e. 400 GAmps) prior to any reconnection. This is 
confined to a current sheet whose width at the coronal point (B' ± = 34GMm~ 1 ) is A ~ 14 
Mm. It is natural that the half-width of this vertical sheet is close to the height of the point 
since the current sheet must totally separate the new and old fluxes. Had the absence of 
reconnection across this one separator been the only constraint on the coronal field, the field 
would contain free energy 

/ A^ 2 

AW ~ — A^ = , (6) 

2c 8vrL ' v ; 

which is AW = 5 x 10 31 ergs. Although there are numerous other separators imposing other 

constraints, this crude estimate is consistent with the amount of energy ultimately released 

in the flare. 

While the width of the sheet, A, depends on global properties of the field, its thickness, 
S, depends on non-ideal processes significant inside the sheet. For a Sweet-Parker solution in 

— 1/2 

the presence of uniform resistivity, rj, the thickness would be 5 ~ AS A where the current 
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sheet Lundquist number is 



Av A ,± |/| i sp 



77 J sp c V 4vr 

Spitzer resistivity in a coronal plasma (T = 10 6 K and n e = 10 10 cm -3 ) yields a characteristic 
current J S p — 10~ 2 Amps, for which the Lundquist number will exceed 10 12 . In this case the 
equilibrium width of the current sheet would be 5 ~ 14 m. Sweet-Parker reconnection at such 
a high Lundquist number transfers flux very slowly, making it essentially "non-reconnection" . 
In spite of the large current density the energy dissipation, P (I / c)d^f / dt ~ AW/r sp , is 
negligibly small. Thus we would expect signatures of a pre-flare current sheet to be extremely 
weak, if observable at all. This is constant with the absence of any current-sheet signature 
from pre-flare EUV images. 

Given this slow reconnection it is possible for the current sheet to exist for any portion 
of the 50-hour build-up phase. The flare itself must arise from a rapid transfer of flux 
across the current sheet. In other words, the current sheet is necessary but not sufficient for 
reconnection. 

In light of the coincidental appearance of numerous thin coronal loops of the new con- 
nectivity it seems that flux transfer was unsteady and spatially intermittent rather than a 
steady, extended electric field along the separator current sheet. Indeed, it has become rec- 
ognized that localized e lectri c fields are a nec essary condition for fast magnetic reconnection 



Biskamp and Schwarzl ( 1200 If ): iKulsrudl (120011 ) 



4. A Model of Time-Dependent Localized Reconnection 

The localized reconnection will occur across a current sheet where field lines of different 
connectivity come into close enough proximity to be interconnected by localized processes. 
As a model of the process consider a Green-Syrovatskii current sheet with a guide field B g 
such as that from Figure [TTJ Such a current sheet separates field lines of the same magnitude 
but differing in angle by AO = 2ta.n~ l (B±/ B g ). The strength of the reconnecting component 
depends on the width of the current sheet, B± = B' ± A/2, so the reconnection angle is 



A9 = 2 tan 



-1 



A 



2B g /B> ± 



The separator in our pre-flare field is characterized by a dimension B g /B'j_ ~ 10 Mm. When 
the current sheet width reaches A = 14 Mm the reconnecting field angle is AO = 70°. 
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The computation above reveals that current sheet reconnection occurs in the presence of 
a significant guide field. A transient and spatially localized reconnection event in such a cur- 
rent sheet will produce a flux tube as illustrated in Figure IT2"1 While all field lines adjacent to 
the current sheet were initially in equilibrium, the reconnection produces sudden disequilib- 
rium from very sharp bends in the new tubes. These bends pro duce a rapid retraction of th e 
tube; each tube pulls itself straight by s liding through the sheet iLinton and Longcopd ( 120061 ) ; 
Linton. Devore. and Longcopd (120091 ) ; lLongcope. Guidoni. and Lintonl (j2009bl ) . The retrac- 
tion is governed by dynamics independent of the process which created the non-equilibrium 
tubes. 




PSP 



v v. 

re 9 ion 1 region 2 




w(tj 



Fig. 12. — Patchy, transient reconnection in a current sheet. Reconnection within patch R 
creates two flux tubes (light and dark grey) shown in the left panel. Each of these retracts 
through the sheet; the central panel illustrates the retraction of the lower tube. Rotational 
discontinuities (RDs) propagate along the tube at the Alfven speed. Central compression 
creates two gas dynamic shocks (GDSs) confining a hot, dense plug of post-shock plasma 
(PSP). This moves downward with the flux tube at speed v z and expands horizontally at the 
shock speed v s . The right panel shows the plug at successive times culminating in t — th, 
when it reaches the Y-type neutral line along the bottom of the sheet. 



4.1. Structure of reconnection outflow 



Post-reconnection retraction releases magnetic energy by shortening the tube on the 
Alfven time characteristic of magnetically-driven dynamics. The initial bend decomposes 
into two rotational discontinuities (RDs) propagating at the Alfven-speed and two slow mode 
shocks (SMSs) or gas-dynamic shocks (GDSs). The RDs produce the rapid shortening, but 
they do so without changing temperature or density of the plasma. This rapid shortening 
takes the form of Alfven speed motions vi directed partly inward in order to decrease the 
length. This parallel compression is therefore much faster than the slow magnetosonic waves 
which govern its dynamics, and the flows collide in SMSs or GDSs. 
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Longcope et al. (l2009bl ; hereafter called LGL09) presented a three-dimensional model of 
such transient, patchy magnetic reconnection. They used thin- flux tube dynamics to describe 
the retraction and compression. They found dynamical evolution following the diagram in 
the center panel of Figure [J2], with purely hydrodynamic GDSs formed by the converging 
flows. 

In the model of LGL09 each RD in the downward-moving tube produces a flow, desig- 
nated region 1, directed along the bisector of the bend, 

Vl = ^2v A0 sin 2 (A#/4)x-t; A>0 sin(A#/2)z (9) 

where va,o is the Alfven speed in the background magnetic field outside the current sheet. 
The two RDs direct flows inward in opposing horizontal senses. In a reference frame moving 
downward with the horizontal segment each inflow has gas-dynamic Mach number (LGL09) 



Mi = ™ = J-^sin 2 (A0/4) , (10) 
c s ,i V 7Po 

where c s> i is the sound speed in the pre-reconnection (and pre-shock) plasma. Except in 
cases of very small reconnection angles AO, the small value of pre-reconnection /3, denoted 
Po, ensures very large inflow Mach number. 

The collision of these two inward flows brings the plasma to rest at two gas-dynamic 
shocks (GDSs). The density increase across these shocks follows from the gas-dynamic 
conservation laws (LGL09) 

p 2 y /9 + AMf + 2M 1 

r = — = = , (11) 

Pi y/§ + AMI ~ M i 

assuming the shocks are in steady state (we have used 7 = 5/3). The temperature ratio 
across the shock 



T2 = 1 
T\ r 



1 + |Mi (Ji + 4mJ/9 + 2Mx/3^ 



12) 



depends on (3 and A9 through eqs. (fTDl and (fill) ; this dependence is plotted in Figure [T3l 



4.2. Alternative reconnection models 



In the discussion above the heating follows from non-resistive, post-reconnection dy- 
namics alone. It is therefore natural that the model of LGL09 gives results similar to other 
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Fig. 13. — Properties of the PSP vs. the reconnection angle A6 fo r various models. Curves 
show n are solid: Longcope et al. (LGL09), dashed: Soward ( 119821 ) and diamonds: Lin and 
Lee (jl994j ). The lower panel shows the temperature ratio, T 2 /7i, for various values of the 
initial plasma-/?, noted on the right. All models agree for small angles for each value of /3q. 
A triangle denotes the situation for a current sheet with A6 = 70° and /3q = 0.002. The 
top panel shows the fraction of the tube between RDs occupied by the PSP for the case of 
(3 = 0.01 for all three models. 
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model s of fast magnetic reconnection. The most well-known such model, that of Petschek 



( 11964 ). is of two-dimensional, steady-state reconnection between perfectly anti-parallel field: 



AO = 180°. In that mod el two long steady outf l ow jet s are confined between standing SMSs 



It was realized early on iPetschek and Thornei (119671 ) that the addition of a guide field de- 



manded a more complicated shock structure involving both RDs and SMSs. For /3 « 1, 
or Ad ~ 180°, the shocks are close together and the configuration appears similar to the 
classic Petschek one. In other cases, such as our own with <C 1 and a significant departure 
from A8 = 180°, the SMSs occur well inside the bends in the field: the RDs. From the 
conventional end-on perspective the configuration still resembles the Petschek model, but 
from the side-view perspective any given field line looks like the flux tube in the center panel 
of Figure [12 but with SMSs instead of GDSs. 

A complete quantitative tr eatme nt of Petschek reconnection in the presence of a guide 



field w as presented by Soward (119821) a nd then re-derived in simpler form by Forbes et al. 



( 119891 ) and Vrsnak and Skender (120051 ) . This analysis reveals that the post-reconnection 
temperature depends on angle and Alfven speed alone, in a manner very similar to the 
transient three-dimensional LGL09 model (see Figure [T3"j) . It does not, for example, depend 
on the reconnection rate or even the microphysics generating the reconnection electric field. 
The distance along a given field line between the SMSs is some fraction of the distance 
between the RDs. In the steady-state, two-dimensional model this fraction depends on A9 
in exactly the same way as the post-shock plasma (PSP) size does in the transient three- 
dimensional model (see top panel of Figure [T3"j) . 

Properties of the more extensively studied Petschek model can be read off the far right 
side of Figure [TBI AO = 180°. In that case the two shocks, RD and SMS, combine into the 
single shock with no gap — the PSP occupies 100% of the region, where the dashed curve 
intersects the right axis of the upper panel. The final temperature ratio depends on f3 ; the 
oft-used value (3 = 0.01 produces a 40-fold increase in temperature, in agreement with the 
dashed curve in the lower panel. 

Figure [13] reveals remarkably good quantitative agreement, especially for for A9 <; 
120°, between seemingly different reconnection models. The level of disagreement at large 
angles results from significant reduction in field strength across the SMSs assumed absent 
from the GDSs of LGL09. This reduction drives the post-shock /3 to non-negligible values 
contradicting the assumption of LGL09 that ^ < 1 on both sides. These assumptions 
do, however, appear to be good ones for smaller angles. In particular, the SMSs do not 
significantly decrease the magnetic field strength when A9 ^ 120°, in contrast to their 
switch-off character in the traditional Petschek limit (Ad = 180°). 



The agreement between a steady two-dimensional model and a non-steady three-dimensional 
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model can be understood in terms of a non-steady one-dimensional model. Lin and Lee 



( 119941 ) solved the Riemann problem beginning with an infinite current sheet bending an 
otherwise uniform magnetic field. This sheet decomposes into a pair of fast magnetosonic 
rarefaction waves, outside a pair of RDs, outside a pair of SMSs. It has been demonstrated 
that this solution is similar to Petschek reconnection in two-dimensions either with or without 



a guide field iLin and Led ( 119991 ) . The diamonds in Figure [T3l show that the SMS separa- 



tion ( relative to the RD separation) and the post shock temperature match those of Soward 



( 119821 ) over all angles. Even a steady, magnetized flow contains individual flux tubes each 
undergoing time-dependent evolution. This evolution obeys non-resistive MHD, independent 
of the reconnection mechanism. Evidently it is equally well described in one, two or three 
dimensions. 

We conclude that post-shock temperatures from fast reconnection across a current sheet 
depends on the field strength and angle according to a fairly robust relationship given by eq. 
(Fl2|) and shown in Figure [131 The specific details of the reconnection appear to be irrelevant 
in determining this particular property. The reconnection outflow may take the form of a 
long jet, as many models have assumed, or it may consist of numerous isolated plugs. Either 
situation produce very similar post-shock properties. For the initial angle A9 = 70° to heat a 
plasma from active region temperature, To = 3 MK, to Ti = 37 MK observed in the loop-top 
source, requires reconnection at /3q ~ 0.002. In a field strength of B = B g j cos(35°) = 400 
G, a 3 MK has this value of beta for n c = 1.5 x 10 10 cm -3 . 



5. A Model of the Loop-top Source 



We suggest that the loop-top thermal emission observed by RHESSI is from a post-shock 
plasma of the kind just described. The numerous distinct loops observed in 17lA suggest 
a spatially intermittent (patchy) reconnection has formed isolated high-temperature plugs 
of the kind shown in Figure [121 The flux tubes would be formed by reconnection across a 
current sheet like that in Figure [TD1 Each plug has the same cross-sectional area as the tube 
and extends a distance w(t) along its axis, as in Figure [121 With flux 5<p and field strength, 
B, set by the flux layers between which it is confined, the flux tube has cross sectional 
area 8ip/B. This is independent of the actual thickness of the current sheet across which 
the reconnection occurs; in f act the LGL09 model and the simula tion which motivated it 
Linton and Longcopd (120061 ); iLinton. Devore. and Longcopd (120091 ) take the sheet to much 
thinner than the tube. The emission measure of this one plug is 



SEM = n 2 e>2 w8ip/B , 



(13) 
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where n e ^ is the post-shock electron density. Provided the plasma is sufficiently dense to 
thermalize, all of the plug's emission is characterized by the post-shock temperature T 2 . 

New patches of reconnection, distributed over the current sheet, will produce new flux 
tubes, two per patch, at an average rate 2<&/5<p. Each of these includes a central plug which 
emits for a time r p i ug . Emitting plugs have a combined emission measure 

EM = 2< 2 ( W )r plug ^ . (14) 

-DO 

where (w) is the axial extent averaged over the emitting life-time. The emission measure 
therefore depends, through properties on the right hand side, on the angle A8 between 
reconnecting field lines and the mean reconnection rate. This rate, $ is not related to the 
reconnection electric field within a single patch, but is rather the rate at which new patches 
are produced throughout the current sheet. 



5.1. The size and life of the PSP plug 

The post-shock plasma (PSP) plug is confined by supersonic inflows generated by the 
rotational discontinuities (RDs) shown in Figure H2J The size, (w) and life, r p i ug , of a single 
emitting plug can be estimated by considering how long these flows are able to confine it. 
This follows from a detailed analysis of the plasma dynamics inside the retracting flux tube 
shown in Figure [121 

We hereafter assume reconnection at sufficiently low /3q (pre-reconnection) and at suffi- 
ciently large angle, AO, that the flow in region 1 is hypersonic: Mi ^> 1 according to eq. (fTUj) . 
Using this in eq. (II ip gives density ratio near the gas-dynamic limit r ~ 4. We hereafter use 
this value to simplify the analysis; using the full expression would produce more complex 
expressions but similar values. The accretion of the inward flows cause the PSP to grow as 
the GDSs propagate at 

v a = ^ ^v xA . (15) 

r — 1 6 

Vertical motion of the tube segment persists at roughly constant speed, v Zj i = —va,o sin(A6 l /2), 
until it reaches the bottom of the current sheet. This happens at t h = h/\v Z) x\ (see Figure 
fT2i) . at which time the PSP occupies 

w(t h ) = 2v s t h = §/itan(A0/4) , (16) 

of the total segment length, 2/icot(A6 l /2). This is the fraction whose dependence on AO is 
plotted in Figure [131 In purely two-dimensional models, i.e. A6 = 180°, the retraction is 
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halted at a fast magnetosonic termination shock Forbes and Priestl (119831 ); here it is a single 
event in the evolution of the tube and the PSP. 

The end of the RDs means the end of the supersonic inflows confining the PSP. The 
vertical forces which halt the retraction do not, however, directly stop the horizontal flow 
throughout region 1. Instead, inflow stoppage propagates horizontally as a rarefaction wave 
(RWi in Figure IT4l) moving at Vf = c Sj i + v x ,\ ~ v Xj i. At time th the rarefaction wave and 
the shock are separated by a horizontal distance 



Ax = hcot(A0/2) - ~w{t h ) = %at(A0/2)-itan(A0/4)] 



(17) 



as shown in Figure HH They move toward one another at relative speed, Vf + v s , meeting 
after 

At f = AX = -^-[3cot(A0/2) -tan(A0/4)l . (18) 
v XtX + v s Av xA 

At this time the shock has propagated a horizontal distance 

h 



UK, 



2v s (t h + At 



f) 



2sin(A0/2) 

since the reconnection episode; this is the maximum extent of the plug (see Figure [1 

Ax ^ 



(19) 
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Fig. 14. — A schematic of the time evolution of the post-shock plasma plug depicted with 
time running downward. The three graphs on the right show density vs. x at characteristic 
times. Regions 1 and 2 are light and dark grey regions confined by shocks and rarefaction 
waves moving at different horizontal velocities. The panel on the left shows the evolution of 
these regions and the different shocks. Horizontal arrows show the extent of the PSP plug 
at specific times. 



The interaction of the rarefaction wave and the gas-dynamic shock results in a new 
rarefaction wave (RW2) propagating toward the center of the plug at the post-shock sound 
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speed c S) 2- This represents the free expansion of the unconfined hot plasma which ultimately 
disassembles the plug. These second rarefaction waves will meet after an interval 



2c 



a, 2 




v x ,i 2^5 



(20) 



marking the end of the PSP. The final expression in eq. f )20|) uses the high Mach-number 
limit of eq. (JTJJ, T 2 /T l (5/9) Aff. The total lifetime of the plug 



plug 



2 '"max 



1 1 



3w„ 



1 ' 

V5. 



(21) 



4tu, sin 2 (A#/4) 

is longer than its free expansion by a factor c Si2 /v s ~ v5. Since the plug expands linearly 
and then contracts linearly its average extent is one-half its maximum: 

h 



(w) 



2 '"max 



4sin(A#/2) 



(22) 



In the analysis above the life time of the emitting plasma was fixed by the time-scale 
for pressure confinement. Several previous inv e stigat ions have used instead a time-scale for 
dissipation by conductive cooling Ijiang et al\ (120061 ) . A detailed analysis by Guidoni and 



Longcope (2010) shows that conductivity produces extended thermal fronts in advance of the 
GDSs. Temperature increases gradually over the fronts reachi ng the steady-state Rankine- 
Hugoniot value at the density jum p which is the actual GDS ( lKennellll987t Ku and Forbes 
19921 ; iGuidoni and Longcopell2010L more accurately called an iso-thermal sub-shock). There 
is very little density enhancement within the fronts, so they would not contribute appreciably 
to the emission measure. Thermal conduction could precede the rarefaction wave RW2, 
causing At c to be even smaller than our estimate in eq. (1201) . This is the smallest contribution 
to Tpi U g so we expect that eq. (12T1) . based on pressure confinement alone, is probably close to 
the actual life. 



5.2. The loop-top emission 

Assuming the reconnection occurs near the center line of the current sheet, of width A, 
the net retraction distance is h = A/2. The current sheet width is related to the reconnection 
angle by eq. (jSJ), 

h = ^tan(A#/2) , (23) 

from the structure of the Green-Syrovatskii sheet with a guide field. The field properties at 
the separator point in Figure [TO] and an angle AO = 70° inferred from RHESSI temperatures 
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gives an average plug size of (w) = 3.0 Mm. Taking the external density to be n e> o = 
1.5 x 10 10 cm~ 3 gives a life-time of r p i ug = 8.5 sec. The loop-top source is somewhat larger 
and considerably longer-lived than a single PSP so it must consist of many such events at 
any one time. 

The net emission measure from the assembly of plugs is given by eq. (fl4j) . Expressing 
the total magnetic field in terms of the guide field, B = B g / cos(A6'/2), and using eq. (|23|) . 
gives 

5/ 2 / i 

EM = 15.4 C '°Y 9 —^7 tt , 24 

B'l sin 2 (A#/4) ' V ; 

where m p = p/n e is the proton mass. In this way the emission measure in super-hot thermal 
material is directly proportional the rate of reconnection. 

There are several ways by which we might derive the reconnection rate $. Equations 
()3]), (JSJ) and OH]), concerning the Green-Syovatskii current sheet, can be combined to yield a 
relation between flux transfer and the decrease in reconnection angle 

. _^A*| _ Bl sin(Afl/2) d(A6) 

dt B' ± cos 3 (A#/2) dt ' 1 ' 

The reconnection angle AO decreases as the current decreases due to reconnective flux trans- 
fer across the sheet. The temperature of the PSP also depends on AO as shown in Figure 
[T3"l It is the decrease in this angle, due to reconnection, which leads to the decreasing tem- 
perature in the super- hot thermal component shown in Figure [9j We can combine these two 
relations to derive A8(t) from T(t) for B g = 331G and an assumed value of pre-reconnection 
density, n e ,o- 

Assuming larger densities n 6i o = 3 - 8 x 10 10 results in slightly larger reconnection angles 
AO = 80° - 100° due to larger fto- The decreasing temperature causes the angle to decrease 
at —Ad/dt ~ 65 degrees per hour for the two density choices. Using these rates and angles 
in eq. f )25]) leads to $ ~ 20 - 50 Gigavolts as shown in Figure [T5j This flux transfer rate 
in eq. f[2"4"]) gives an emission measure commensurate with RHESSI observations when the 
larger initial density, n c $ = 8 x 10 10 cm -3 is used. 

A more direct, and more conventional, measurement of $ can be obtained by map- 
ping the apparent motio n of flare ribbons across th e magnetic field in the lower atmosphere 
Forbes and Priestl ( 19841 ) ; Poletto and Kopp ( 1986). This is frequently done with high ca- 



dence TRACE UV (1600 A) images ISaba. Gaeng. and Tarbelll fl2006h : lLongcope et all (120071 ) 



but these were unavailable at this time since TRACE was imaging instead in EUV 171 A. 
Observations at this wavelength usually reflect >1 MK plasmas in the corona, however, the 
bandpass does include contribution by enhanced chromospheric emission at the feet of flare 
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2004-Feb-26 

Fig. 15. — Plots of inferred flux transfer rates (bottom) and the emission measure they would 
produce. Blue curves are derived from the temperature change using eq. ( 1251) and red curves 
are from the motion of 171 A ribbon features across the magnetogram. Curves with triangles 
(squares) are derived using n e = 3 x 10 10 cm -3 (n e = 8 x 10 10 cm~ 3 ). The black curves 
are the super-hot emission measure derived from RHESSI, error bars are the values from 
detectors 3 and 4 alone. 
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loops especially in the early phase of the flare. Such emission is evident in the 1:59:14 and 
2:03:14 panels of Figure HJ 

We therefore use 17lA images to compute the reconnection flux transfer rate. We 
outline the chromospheric ribbons seen at the start of the flare, and, as the flare evolves, 
track only the outward spread of the ribbon, i.e. the spread of the brightening away from 
the polarity inversion line. This approach avoids contamination due to emission by multi- 
megakelvin plasmas in the post-flare loops, which have cooled down to emit in 171 A between 
the expanding ribbons. 

The flare ribbons are tracked in co-aligned TRACE 171 A images from 1:00 UT for two 
hours with a cadence of approximately 40 s and a pixel scale of 0.5 ". The chromospheric 
ribbons are then mapped in a magnetogram co-aligned with the flare images and the total 
magnetic flux covered by the ribbons is measured. To approximate the magnetic field in 
the chromosphere, we extrapolate the photospheric line-of-sight magnetic field, obtained by 
MDI/SoHO 20 minutes before the flare onset, into a nominal chromosphere height of 2000 
km using a potential field assumption. Num erous expe r iment s have shown that measure- 



ment of the reconnection flux is quite robust iQiu et all (120071 ) . The extrapolation merely 



decreases the measured reconnection flux by 25%, but does not modify the time profile of $. 
Assumptions other than a potential field assumption (such as a linear-force-free assumption) 
and/or extrapolation to a different height < 2000 km would not change the $ profile, and 
its magnitude would be modified by <25%. 

We also estimate uncertainties in the flux measurement by mis-aligning MDI magne- 
togram and TRACE images by up to 2" and by using a set of brightness threshold values 
for ribbon detection. The upper-limit of the uncertainty due to these sources is 30%. Given 
the geometry of this flare, the flare ribbon in the negative magnetic field is more accurately 
tracked. Therefore $ measured in the negative field is used and plotted in red on Figure [T^J 
Using this rate in eq. ( 125]) gives independent measures of the super-hot emission measure. 
Once again we find agreement when the larger initial density, n e ^ = 8 x 10 10 cm" 3 is used. 

Relation ( 123]) could be inverted to use the observed super-hot emission measure to 
estimate $. Agreement in Figure [TBI shows that doing so yields a third independent measure 
in agreement with the other two - temperature decrease and ribbon motion. It appears that 
the flare is driven by flux transfer at a mean rate $«4x 10 18 Mx/s (40 Gigavolts). We 
stress once more that this is not a local electric field, but rather the net effect of multiple 
reconnection events. 
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6. A Model of the flare 

6.1. The loop-top source 

The TRACE loops appearing in the 17lA images (i.e Figure H]) have typical diameters 
of 2 Mm. Attributing to them the same field strength as the separator, B w 330 G, means 
each has flux 5ip = 10 19 Mx. Patchy reconnection at the inferred rate, $ = 4 x 10 18 Mx sec -1 , 
would produce such flux tubes at the rate of 2$/<5<^ = 1 per second. The PSP plug at the 
center of each tube would persist for r p i ug ~ 8 sec, so there would be ~ 8 such plugs visible 
at any time. Using the A8 = 100° angle inferred from the emission measure computations 
means the plugs are confined between |B| = 430 G field, and are thus 1.7 Mm in diameter; 
they range in length from zero to w max = 5 Mm. We posit that the 9 x 16 Mm loop-top 
source in the 25 - 50 keV RHESSI images (Figure [Sb - c) is such a superposition of plugs. 
While 8 plugs are visible at any instant, each image shows a super-position of ~ 28 plugs 
visible for some portion of its 20-second integration (see Figure [TBI) . 

The volume of a super-hot plug would average (w)5<j)/B ~ 5 x 10 24 cm 3 , so the combined 
volumes of the 8 plugs simultaneously visible at one instant is V = 4 x 10 25 cm 3 . Assigning to 
the observed loop-top source a volume V s = 7r(4.5 Mm) 2 (16 Mm) = 10 27 cm 3 , as estimated by 
the 70% contour of 25 - 50 keV images in Figure [T6l it would have a filling factor of / = 0.04. 
(Had we used the 50% contour instead, the filling factor would have been / = 0.01.) A post- 
shock density of n e 2 = 4n 6) i = 3 x 10 11 cm -3 gives the ensemble of plugs an emission measure 
of EM = 4 x 10 48 cm -3 , in good agreement with Figures 191 or [TBI 

Transferring the entire domain flux, A^f^-i = 2.8 x 10 21 Mx, would take r rx ~ 12 minutes 
at a steady reconnection rate $ = 4 x 10 18 Mx sec -1 inferred above. Each reconnection 
event transfers flux 5<p = 10 19 creating two post-reconnection loops — one in each post- 
reconnection domain. Thus we expect ~ 560 flux tubes will be created over w 700 seconds 
(12 minutes). It is this transfer rate, rather than cooling following impulsive heating, that 
sets the overall flare time scale. 



6.2. Post-flare loops 

We have presented a model for how fast reconnection could produce super-hot plasma 
just outside the diffusion region. One plug of this super- hot plasma would persist for r p i ug ~ 8 
seconds before the confining flows cease and its own high pressure caused its disassembly. 
Thermal conduction during and following the plug's compression would eventually reach the 
footpoints of the flux tube causing evaporation. Details of this phase are beyond the scope 



1 :S8:00-1 :58:20 
12-25 keV 25 




Fig. 16. — Images of the loop-top source, (top) A super-position of 28 plugs 2 Mm across 
with lengths varying randomly to 6 Mm. They are distributed randomly in space near the 
separator (magenta curve). Blue lines show the extent of the RD along the field lines, and 
green diamonds show the locations of the GDSs. A red asterisk marks the centoid of the 
super-positon. (bottom) RHESSI images from 12 - 25 kev (inverse grey scale) and 25-50 
keV (green contours) from an integral over 1:58:00 - 1:58:20, at the peak of the thermal 
phase. The red curves show an ellipse with moments matching the central (70%) contour 
of the 25-50 keV emission. The major and minor axes, 16.0 Mm and 9.1 Mm, are straight 
lines. 
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of our model, however, it is possible to predict its overall contribution. 



We suggest that the evaporation from the thermal conduction fills the loop with lower- 
temperature plasma observed as the cooler (T ~ 20 MK) component of the RHESSI spectra. 
Following Antiochos and Sturrock (119781 ) we assume the conductive cooling transfers all 
energy from the thermal and kinetic energy of the tube to the evaporated component. The 
temperature and emission measure of this component matches the plasma observed in X- 
rays by GOES (see Figure Ml ■ Such observations are c ommonly explained in terms of the 
Neupert effect iNeupertl (119681 ) ; iDennis and Zarrol (119931 ) whereby direct flare energy release 
evaporates chromospheric material. We propose here that, at least in this present case, the 
evaporation is driven by thermal conduction rather than non-thermal particles. 

Images from RHESSI's 6-12 keV often resemble other soft-ray images and can be 
used to estimate the extent of the plasma detected by GOES. The 50% red contours Figures 
[H^l d outlines a 30 x 16 Mm region, suggesting an ellipsoid of volume V = 4 x 10 27 cm 3 . 
Using a filling factor, / = 0.18, justified below, we estimate the density of evaporated 
material as n P = y/EM/fV. Using EM from ratio of GOES channels gives a maximum 



n P 



2.6 x 10 11 cm-' 3 at 2:04:12, thereafter falling steadily to n e = 6 x 10 10 at 3:00:00. 



Once a fully retracted tube is full of T « 20 MK plasma from evaporation it will cool 
throu gh the rmal conduction and then radiation. According to the cooling model of Cargill 
et al. ( 119951 ) . a loop at the temperature and density inferred from GOES, and length given 
by eq. (jTJ, would cool to half its temperature in between 1 minute (at 1:58) and 3 minutes 
(at 2:30). They would reach T ~ 1 MK shortly thereafter and appear in the TRACE image. 
This explains loops visible by TRACE almost immediately, even as the temperature inferred 
by GOES approaches 1 MK much more gradually. The continued appearance of the these 
loops over one hour is due to the continued creation of new loops through reconnection. 



6.3. Energetics 

Each tube of reconnected flux 5(j) = 10 19 Mx releases energy by decreasing its length by 
a distance A£ = 2/itan(A0/4). This decrease in length, with no appreciable change in field 
strength, will decrease the flux tube's magnetic energy by 

SE = ~ 10 29 erg , (26) 

8n 

using values inferred for the current sheet, with h = 7 Mm and A8 = 100°, B = 430 G. The 
expansion of the external flux layers into the volume vacated by the retraction will release 
an equal amount of energy (LGL09) leading to a net release of SE mag = 2 x 10 29 ergs. This is 
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the value at the onset of the flare (1:58:00); later flux tubes will release steadily diminishing 
energies as the current sheet shrinks causing A9 and h to decrease. If the average energy 
release is half the maximum value, the complete set of 560 flux tubes would release 5.6 x 10 31 
ergs, consistent will the initial free energy, AW, from eq. (J6j). 

At a time th = 3 sec following its reconnection a given tube has fully retracted and 
released all its magnetic energy. At this instant the plug is with) = 2.2 Mm long and 
contains thermal energy 

U^ lug \t h ) = 3k h T 2 n e , 2 (5<j ) /B)w(t h ) ~ 2 x 10 28 ergs (27) 

only 10% of the magnetic energy liberated by the reconnection, SE mag . This is rather low 
efficiency compared to the anti-paralle l case (A6 = 180° ) wher e 40% of the free energy is 
directly thermalized by the slow shocks iPriest and Forbes! (120001 ) . In our case the other 90% 



of the liberated energy was converted into bulk kinetic energy of retraction. Roughly one- 
sixth (sin 2 [A#/4]) of that remaining kinetic energy is directed horizontally and is ultimately 
thermalized before the plug disassembles itself. Thus 25% of all magnetic energy liberated 
by reconnection appears as thermal energy in the loop-top source at some point. Since a 
given plug will have, on average, half of this amount at a random time, the eight plugs 
simultaneously visible will contain 

= 8x3k h T 2 n ei2 (5(j ) /B)(w) ~ 10 29 erg . (28) 

This is our estimate of the thermal energy contained in the super-hot loop-top source at its 
peak. 

The instantaneous thermal energy content of the low-temperature component, filled 
with evaporated material, is U^ OES ^ = 3k\ ) TyJEM V f , where we take T and EM from 
the GOES channels and assume a filling factor / = 0.18 for reasons described below. The 
thermal energy, shown in Figure [T7J rises to a peak of 2 x 10 30 ergs at 2:01:25. At 1:58:00, 
when we found the super-hot component energy, the low-temperature component has an 
order of magnitude more (red triangle). 

The low-temperature component cools by radiation and by thermal conduction. Ra- 
diative losses, computed using the temperatu re dependent Mewe radiative loss function 



Mewe. Gronenschild. and van den Oordl (119851 ) , are plotted in blue on Figure [TTJ The rate 
of conductive energy loss to the chromosphere is P CO nd = ^th/Veond where 

— ! ! ~ (29) 

r cond 3k h Tn e (L/2) 3A; b n e (L/2) 2 ' 1 ; 

for a loop of full length L; Spitzer conductivity is k = kqT 5 ^ 2 . The EUV loops traced in Figure 
[5] had apparent (plane-of-the-sky) lengths increasing from 7 Mm to 12 Mm according to 
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Fig. 17. — Time histories of the energy (top) and power (bottom) from the flare. Green 
and blue curves are losses from conductivity and radiation respectively. These losses are 
integrated backward to produce cumulative energies in the top panel — the net energy lost 
after that time. The red line in the top panel is the thermal energy content, from eq. 
( )28l) . in the lower panel it is —dll^/dt. The instantaneous thermal energy content of the 
low-temperature and high-temperature components at 1:58:00 are shown in the top panel by 
o and x respectively. 
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empirical relation eq. ([Q). Doubling this time-dependent function to account, approximately, 
for the third dimensions and portions not observed at 171 A gives conductive cooling times 
ranging from r cond = 10 sec at 2:00, to r cond = 280 sec at 3:00. The resulting conductive 
losses are plotted in Figure [T71 These losses actually represent a transfer of energy to the 
lower temperature atmosphere, where it is presumably radiated. This radiative loss is known 
to far exc eed the radiation fr om optically thin coronal plasma inferred from the GOES X- 



ray bands lEmslie et all (120051 ) . We use the conductive losses to approximately quantify this 



other loss mechanism. 

At T ~ 20 MK conductive losses dominate radiative losses by up to an order of magni- 
tude. We chose filling factor / = 0.18 in order that the two beco me equal at late times, desig- 



nated RTV in Figure [13 as they should be in a static equilibrium iRosner. Tucker, and Vaiana 



( 119781 ). Prior to that, s hort conductive cooling times, 10 - 280 se c onds, demand a sequence 



of new, dynamic loops IWarren. Winebarger. and Hamilton! ( 120021 ): IWarrenl (120061 ) . 



Integrating each of the losses gives the cumulative losses shown in the top panel of 
Figure [T7J The total integrals of radiation and conduction, 2.3 x 10 30 ergs and 3.3 x 10 31 
ergs, are shown by a square and diamond, respectively. We estimate the total energy lost 
by the low energy component as the sum of these two AW = 3.8 x 10 31 (asterisk). This is a 
slightly smaller than the MCC estimate of free magnetic energy, eq. (jSJ), suggesting, among 
other possibilities, incomplete magnetic relaxation. 

To ascertain the degree of reconnection we integrate the ribbon-swept flux whose change, 
$ is plotted Figure [15j The flux discrepancy 

/oo 
${t')dt' , (30) 

is the same as the ribbon flux, but accumulated from the final time. The initial value of 
this flux, A® = 3.3 x 10 21 Mx, is 20% larger than A\l>4_^. The later is only a single post- 
reconnection domain, albiet the largest one, so it is perhaps unsurprising that it falls short of 
the total flux swept by the flare ribbons. Using the ribbon-swept flux in the magnetic energy 
estimate, eq. (JBJ), with the same separator L = 54 Mm, gives an estimate of the magnetic 
energy liberated by reconnection. The total, AW = 8.5 x 10 31 ergs, half-again as large as 
the estimate using A^^_i and twice the energy released by the flare. Indeed, throughout 
the thermal phase of the flare, after 1:56:40, the free magnetic energy tracks the total flare 
energy very well. It seems that energy release through magnetic reconnection is a viable 
candidate for producing the observed flare emission. 
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7. Discussion 

We have presented here a new model for the shock heating of plasma as magnetic energy 
is released following magnetic reconnection. The flux transfer achieved by the reconnection 
electric field triggers this energy release but does not directly participate in it. As a result 
its details do not enter the prediction, and we have bypassed their consideration. The 
model presented here applies a novel treatment of transient, localized reconnection in three 
dimensions, but its shock heating is closely related to the traditional two-dimensional, steady- 
state model of Petschek. The role of shocks are similar in both cases, but the former does 
not occur in an out-flow jet geometry. 

We have applied this model to RHESSI, TRACE, MDI and GOES observations of a 
particular compact X-class flare. Magnetic energy is stored during a 50-hour emergence of 
flux within the active region, and then released by reconnection over less than one hour. 
The majority of energy released in the flare appears to occur through the production of a 
super-hot (T ^ 30 MK) thermal loop-top source, rather than through the acceleration of 
non-thermal particles. Our model predicts temperature and emission measure of the source, 
roughly consistent with those inferred from fitting the RHESSI spectra. 

The observed source is modeled as a super-position of short-lived plugs of plasma heated 
and compressed by the post-reconnection shocks. While this ensemble never contains more 
than 1% of the magnetic energy ultimately released by the flare, it is the first stage in the 
entire energy release process. The subsequent expansion and conductive cooling of the plugs 
drives chromospheric evaporation to produce the more plentiful low-temperature plasma 
(T <^ 20 MK) observed in soft X-rays. These cool freely to appear minutes later, when 
T ~ 1 MK, in TRACE EUV images. The numerous discrete loops are a direct consequence 
of the transient, localized reconnection which creates them; their early appearance is a result 
of the rapid, independent evolution of each loop. The smooth evolution of the ensemble is a 
result of the gradual decrease in current as flux is transferred. 

We assumed the reconnection occurred within a current sheet which had formed at some 
time prior to the flare. Current sheets are pre-requisites for traditional fast reconnection 
models since they bring field lines of significantly different connectivity close enough to be 
mutually affected by a small-scale reconnection electric field. We estimated the properties 
of the current sheet assuming it be in equilibrium or reconnecting slowly according to the 
Sweet-Parker model. (Any possible resistive dissipation in this pre-flare current sheet is 
neglected in our energy-release scenario.) The two current sheet properties most significant 
for our model are the strength of the equilibrium magnetic field adjacent to it and the angle 
between the magnetic field it separates. 
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In addition to the magnetic properties, our model depends on the density of the flux 
tube plasma prior to its reconnection. This density determines both the Alfven speed and 
the post-shock density. We found that large observed emission measures demand fairly large 
densities, n e0 ~ 8 x 10 10 cm~ 3 . This demand is not unique to our model, since it derives 
from the limit on density enhancement at a single shock. Indeed, our oblique shocks, similar 
to gas-dynamic shocks, produce a larger enhancement factor than switch-off shocks of two- 
dimensional models: 4 versus 2.5. It thus seems impossible that the super-hot post-shock 
material could have arisen from pre-shock material at densities comparable to those observed 
in non-flaring active regions, n e w 10 9 . Something in the flare process evidently "pre-fills" 
the flux — at least the flux immediately adjacent to the current sheet. 

While the high post-shock density was inferred from the observed emission measure, 
independent evidence for it comes from the occurrence of thermal bremsstrahlung from 
a plasma at T w 40 MK confined to <^ 10 Mm along a magnetic field. The mean-free 
paths of electrons at that temperature, £ m f p = 100Mm(10 9 /n ei 2), will greatly exceed the size 
of the source when density is low. Low density plasma streams accelerated by field line 
retraction (i.e. the RDs) would presumably pass through one another, rather than shocking. 
Such counte r -strea ming has been observed in simulations of collisionless fast reconnection 



Drake et al\ (120091 ) . The fact that we observe bremsstrahlung from a Maxwellian electron 
population implies sufficient collisionality to thermalize the loop-top source: £ m f P ^ 1 Mm, 
and therefore very high post-shock density; n e ^ ^ 10 11 cm -3 . 

The pre-reconnection density n e0 is not, strictly speaking, a pre-flare density. The 
thermal phase of the flare, to which our modeling has been devoted, began three minutes 
after the flare itself. We have no estimate of the coronal density, even post-reconnection, 
from the first three minutes (1:53:30 — 1:56:40) where thick-target emission was the best fit. 
It is entirely possible that the pre-flare density was far lower than the density we infer for flux 
tubes reconnected during the thermal phase of the flare. Indeed, all observational evidence 
suggests pre-flare densities typical of ARs: n e w 10 9 . One possibility for large n e) o during 
that phase is that energy release during the initial non-thermal phase somehow enhanced 
the density of flux tubes in the vicinity of the current sheet prior to their reconnection. 

We have not attempted to address the non-thermal phase with our purely fluid model, 
and therefore cannot include its energetic contribution in our estimate. It is over this period 
that the inferred magnetic energy input diverges most from the energy output in Figure 
[T71 Our output rate neglects the precipitation of non-thermal particles, and uses conductive 
losses assuming collisional, Spitzer transport. We therefore expect the input and output to 
disagree to some extent during this period, as they do. 



In contemplating possible sources of non-thermal particles it is worth recalling the failure 
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of fluid model when pre-reconnection densities are too low. Low collisionality in the post- 
shock density leads to a non-thermal electron population with mean free paths comparable to, 
or larger than, the total loop length. Furthermore, higher pre-reconnection Alfven speeds 
lead to even higher post-shock temperature predictions, T 2 ~ v\ ~ n~Q. While these 
extremely large temperatures could not be collisionally thermalized within the loop, they are 
based on energy conservation laws which must still obtain. We therefore expect a population 
of high energy, non-thermal electrons, accelerated directly by the Alfvenic loop retraction, 
the details of which awaits future kinetic calculations. 

The above reasoning suggest that the non-thermal phase occurs when the pre-reconnection 
densities are low. The ambient, pre-reconnection densities rise, perhaps due to evaporation 
driven by observed precipitation of non-thermal particles, and subsequent reconnection oc- 
curs at densities high enough to produce collisionally thermalized plasma at the loop-top. 

This evaporation scenario is similar to standard flare models except that ours requires 
evaporation to occur on a flux tube before that flux tube reconnects to release its energy. 
It is the rapid energy release that raises the temperature of the already-dense material 
to T m 40 MK. Evaporation flow is unlikely to cross field lines so evaporative pre-filling 
seems to require a single flux tube to undergo multiple reconnections. In fact, models of 
reconnection in more complex magnetic topologies do predict a single flux t ube to reconnect 



severa l times, passing through dif f erent separators into intermediate domains ILongcope et al 



( 120071 ); ILongcope and Beveridgd (120071 ) . Moreover, we observe post-flare loops (Figure HJ) in 
domains, such as P09 - iVOl, before we observe them in P04 - iVOl, for which we compute 
the net flux transfer. Sufficient time between successive reconnections would permit a tube 
to be "pre-filled" by evaporation driven by energy released by this preceding reconnection. 

The requirement that evaporation pre-fill loops before their final, energy-releasing re- 
connection, offers insight into the rarity of thermal flares of the kind we have modeled. 
Evaporative pre-filling is probably easier, and therefore a more common occurrence, in small 
compact flares, with short loops, than in large eruptive flares. The loops in this flare, ranging 
from L = 13-22 Mm (doubling the visible EUV lengths in eq. p]), are short compared 
to those in eruptive events. Furthermore, the reconnection was estimated to occur very low 
in the corona: « 10 Mm. The emergence of new flux into this small active region appears 
to have loaded it with AW ~ 10 32 ergs, which was later liberated by reconnection. This 
produced short post-reconnection loops which were quickly filled even as they continued to 
be shortened. 

Uncommon are flares whose energy release is sufficiently dominated by thermal emis- 
sion to be modeled only as a fluid, however, we believe our energy release scenario may 
have broader applicability. Large eruptive flares often exhibit supra-arcade downflows, 



-42 - 



which offer evidence for the s ame kind of post-reconnection ret r action used in our mo del 



McKenzie and Hudson! (119991 ) ; ISheeley. Warren, and Wangl ( )2004j ); lMcKenzie and Savage! (120091 ) 
Ultimately any kind of magnetic reconnection releases energy by permitting field lines to be- 
come shorter (except in the special case of anti-parallel fields, where some of the field energy 
is annihilated). If it is fast reconnection then the shortening will occur at the Alfven speed, 
and plasma on the field lines will be compressed super-sonically. We therefore expect the 
energy release scenario in our model to be fairly common, although perhaps not always treat- 
able by fluid models when densities are too low. Even in large, two-ribbon flares, however, 
the gradual phase shows little evid ence of non-thermal particles and yet can account for sub- 
stantial share of the energy release lEmslie et al\ (120051 ). Since it persists far longer than the 
cooling time of a single loop, the gradual phase must be maintained as a succession of new 
reconnection sites create a nd energize new loops IWarren. Winebarger. and Hamilton! (120021 ) ; 
Reeves and Warren! (120021 ). The energy release in these cases may well occur through the 
same shock heating scenario introduced here, evidently without any non-thermal particles. 
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